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Abstract. Recent results have questioned the description of the QSO 
luminosity function in terms of a pure luminosity evolution and call for 
a luminosity dependent luminosity evolution. Measurements of the QSO 
clustering amplitude and evolution allow further distinguishing among the 
various physical scenarios proposed to interpret the QSO phenomenon. 
The general properties of the QSO population would arise naturally if 
quasars are short-lived events connected to a characteristic halo mass 
~ 5 • 10 12 M.q. This is the typical mass of groups of galaxies in which the 
interactions triggering the QSO activity preferentially take place. 



1. Introduction 

QSOs have played an important role as cosmic probes of the young universe: they 
are used as light beacons for absorption-line and gravitational lensing studies, 
as markers of galaxy formation activity and have guided astronomers in hunting 
up primeval galaxies. Their very existence at high redshifts is a challenge for 
cosmological models and historically the evolution of the QSO population has 
been one of the first evidences for an evolving universe. The nature of QSOs and 
the causes of their evolution, however, are far from being fully understood. In 
the standard "demographic" approach the basic physical ideas are tested against 
the observed shape and evolution of the QSO luminosity function (LF). As usual, 
with the growth and improvement of the databases this type of interpretation 
has become significantly more complex and we will argue that to disentangle 
the physical evolutionary patterns of QSOs additional information is needed, 
for example about their clustering properties. 



2. The QSO Luminosity Function 

In the early nineties the general consensus among the astronomers was that the 
shape of the QSO LF is a double power-law and its evolution, at least up to 
z ~ 2.2, can be satisfactorily represented by a pure luminosity evolution (PLE) 
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Figure 1. The QSO number counts for 0.3 < z < 2.2. The continu- 
ous line represents the best fit LDLE model, the dashed line represents 
our best fit PLE model, while the dotted line represents the best fit 
PLE model by Boyle et al. (1992). 



(Boyle et al. 1988): 

$(M B ,z) = 10 o. 4 [ Ms _ Ms ( 2 )]( Q+1 ) - 1 q0.4[M s -M s ( Z )](/3+1) W 

where a and (3 correspond to the bright-end and faint-end slopes of the optical 
LF, respectively. With this parameterization the evolution of the LF is uniquely 
specified by the redshift dependence of the break magnitude, Mb(z) = M* B — 
2.5 k log(l+z), corresponding to a power-law evolution in the optical luminosity, 
L* oc (1 + z) k . Such a result was essentially based on the AAT (Boyle et al. 
1990) and PG (Schmidt & Green 1983) surveys with the addition of the faint 
data in the Marano field (Zitelli et al. 1992). The relevant parameters were 
estimated to be a = -3.9, (3 = -1.5, M* B = -22.4, k L = 3.45 with a redshift 
cutoff z max = 1.9 after which the evolution stops. Then a new generation of 
bright surveys came: the LBQS (Hewett et al. 1995 and references therein), 
HBQS (Cristiani et al. 1995), EQS (Goldschmidt et al. 1992), HES (Wisotzki 
et al. 1996). It turned out that not only the surface density of bright QSOs 
(B ~ 15) is definitely higher than it was thought before, but also the new LF 
is qualitatively different and cannot be reconciled with a PLE (Fig. ^). The 
flattening of the slope in the bright part of the low-z LF is apparent and in the 
redshift bin 0.3 < z < 0.6 the shape becomes consistent with a single power-law 
(Goldschmidt & Miller 1998). This behavior can be modeled with a luminosity 
dependent luminosity evolution (LDLE), in the form of an evolution of the break 
magnitude: M B (z) = M* B (2) - 2.5Hog[(l + z)/3] where M* B {2) is the break 
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Figure 2. The QSO luminosity function. (Left) The dashed 
lines correspond to the best-fitting Pure Luminosity Evolution with 
a = -3.7, f3 = -1.4, M* B = -22.3, k L = 3.26. (Right) Fit with a 
Luminosity Dependent Luminosity Evolution model (see text). The 
points connected with a continuous line represent the observations. 
Error bars are based on Poisson statistic and correspond to 68% (la) 
confidence intervals. 

magnitude at z = 2 and 

forM B <M B (z) : k = k x + k 2 (M B - M B (z))e~ z/A0 
for M B > M B (z) : k = k\. 

In this way the best-fitting parameters turn out to be a = —3.76, (3 = —1.45, 
M* B (2) = -26.3, fci = 3.33, k 2 = 1.17 (La Franca & Cristiani 1997). How robust 
is this result? Could this departure from PLE be due to observational biases? 
The new bright surveys have been carried out with different techniques: UVx, 
slitless, and also cross-checked with other approaches, e.g. the variability method 
(Hawkins & Veron 1995), and the results are essentially matching. Effects due 
to the photometric errors (the Bennet bias), the non-unique SED of QSOs, 
the variability, the representation of the LF in discrete redshift bins are also 
accounted for (La Franca & Cristiani 1997). However, other biases could be 
lurking in: obscuration effects, either intrinsic or intervening (Fall & Pei 1993, 
Webster et al. 1995); effects due to the presence of the host galaxy, especially 
for low-z, low-luminosity objects, although the authors of various surveys state 
that this latter problem has been considered and precautions taken accordingly. 
In this respect, it will not be useless to compare the optical results with radio 
and x-ray data. 

Studies of radio QSOs show that there is an increase of the radio-loud 
fraction among the bright (M B < —25) QSOs, both in x-ray and optically 
selected samples (La Franca et al. 1994, Hooper et al 1996). Possibly there is 
also an increase of the radio-loud fraction at lower redshifts, but this result relies 
entirely on the PG sample and should be considered with caution. There are 
evidences in favor of a slower evolution of the optical LF (OLF) of radio-loud 
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QSOs, with a 2.7 < Url < 3.1. Hook et al. (this conference, Fig.6) show that 
the space density evolution is slower for more powerful radio sources. In this 
way it becomes possible to interpret the flattening of the bright part of the OLF 
as a slower evolution of bright, radio-loud QSOs. 

X-ray QSOs, on the other hand, seem to follow an evolution pretty similar 
to the optically selected ones, with k x ~ 3.2 (Boyle et al. 1994), which is not 
surprising, given the relationship between the optical and X-ray luminosity of 
quasars (La Franca et al. 1995). Recent deep and shallow ROSAT surveys 
provide evidences for departures of the XLF from a PLE pattern (Hasinger 
1998) and there might even be an indication of flattening of the XLF at low-z 
and high luminosities (see Fig. 1 of Hasinger 1998). 

3. Physical Considerations on the Luminosity Function 

At present there is a broad consensus that the QSO phenomenon is the result 
of accretion onto super-massive black holes (BH) in the nuclei of galaxies and 
a number of authors have linked the change in the QSO activity to a corre- 
sponding change in fuel supply at the center of the host galaxies (Cavaliere &: 
Szalay 1986; Wandel 1991; Small & Blandford 1992; Haehnelt & Rees 1993; Yi 
1996). The epoch of QSO formation is related to the time when the first deep 
potential wells form in plausible variants of hierarchical cosmogonies and var- 
ious mechanisms for the evolution of the accretion modes have been explored 
(see Haehnelt, Natarajan and Rees 1998, Cavaliere this conference) to reproduce 
the observed luminosity function and its evolution. For example Cavaliere et al. 
(1997 and this conference) model the rise and fall of the QSO LF as the effect 
of two components: the newly formed BH, which are dominant at z > 3 and 
the reactivated BH, which dominate at z < 3. The reactivation is triggered 
by interactions taking place preferentially in groups of galaxies, of typical halo 
mass ~ 5 • 10 12 M Q . In this scenario the flattening of the bright end of the LF 
at low-z arises naturally due to the re-emergence of the first component which 
is overwhelmed by the second at intermediate redshifts. 

In general, the agreement obtained by the various authors with the observed 
LF within the supermassive BH paradigm and also outside it (e.g. Boyle and 
Terlevich 1997) is remarkable. The study of the QSO clustering provides a 
powerful tool to unravel the riddle of the interpretations. 

4. Clustering of QSOs 

QSO clustering was detected for the first time more than ten years ago (Shaver 
1984). At the beginning the measurement was based on inhomogeneous catalogs 
but soon it became possible to use the same statistically-well-defined samples 
used for the LF works, although with a significantly reduced number of QSOs. 
The two-point correlation function (TPCF, Peebles 1980) has been applied as a 
standard approach to investigate the QSO clustering, but other techniques have 
also been explored, such as the minimal spanning tree (Graham et al. 1995), 
fractal analysis (Andreani et al. 1991), counts in cells (Andreani et al. 1994). 

In recent times complete samples totaling about 2000 QSOs have been used, 
providing a 4 — 5a detection of the clustering on scales of the order of 6/z -1 co- 
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Figure 3. The integral correlation function £(r) - defined as £(r) = 
\ Jo x2 ^{ x )dx - for the quasars in the SGP sample in two redshift 
ranges 0.3 < z < 1.4 (open circles), and 1.4 < z < 2.2 (filled circles). 



moving Mpc (Andreani & Cristiani 1992, Mo & Fang 1993, Croom & Shanks 
1996). The evolution of this clustering is not clear. An amplitude constant in 
comoving coordinates or marginally decreasing with increasing redshift has been 
suggested, an amplitude which appears to be consistent or slightly larger than 
what is observed for present-day galaxies and definitely less than the clustering of 
clusters. Measurements of the quasar-galaxy correlation indicate an amplitude 
about 4 times the galaxy-galaxy clustering (Fisher et al. 1996). With all the 
caveats due to their exotic nature, QSOs display a number of appealing proper- 
ties when compared to galaxies as cosmological probes of the intermediate-linear 
regime of clustering, on distances > 20 Mpc: they have a rather flat n(z), allow 
defining samples which are in practice locally volume limited, their point-like 
images are less prone to the surface-brightness biases typical of galaxies and 
they sparse-sample the environment. The sparse sampling is an advantage, in 
the sense that the contribution of each QSO or QSO pair to the correlation 
functions in practice turns out to be independent, making statistics easy, but is 
also a nuisance, because the sensitivity to the clustering signal is low. 

In an attempt to improve the situation, while waiting for the 2dF QSO 
redshift survey, we have carried out a survey in the South Galactic Pole (SGP) 
over a connected area of 25 square degrees down to Bj = 20.5. Stacked UKSTU 
plates were used to select UVx candidates and the multi-fiber spectrograph 
MEFOS at ESO to take spectra of them. The final sample is made up of 388 
QSOs with 0.3 < z < 2.2 (La Franca et al 1998a). 

The TPCF analysis gives an amplitude r Q = (6.2 ± 1.6) Mpc, in agree- 
ment with previous results. However, when the evolution of the clustering with 
redshift is analyzed, evidence is found for an increase of the clustering with in- 
creasing z (La Franca et al. 1998b). The effect is small, at a 2a level (Fig. ||), 
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Figure 4. The amplitude of the £(15 hr l Mpc) function of z. 
Filled circles: the low- and high-z SGP subsamples; open circle: the 
SGP sample plus the Boyle et al. (1990), La Franca, Cristiani and 
Barbieri (1992), and Zitelli et al. (1992) samples; open triangles: same 
as open circle but divided in two redshift slices; filled triangle: low-z 
AGNs from Boyle and Mo (1993) and Georgantopoulos and Shanks 
(1994); open square: the high-z sample from Kundic 1997. The dotted 
line is the e = —2.5 clustering evolution, the dashed lines are the 10 12 
and 10 13 M Q h' 1 minimum halo masses clustering evolution according 
to the transient model of Matarrese et al. (1997). 



but is interestingly corroborated by other results in the literature. Boyle and 
Mo (1993) measured the clustering of low-z QSOs in the EMSS and Georgan- 
topoulos and Shanks (1994) used the IRAS point source catalog to measure 
the clustering of Seyferts. Altogether a low value of the TPCF at 15 Mpc and 
z=0.05 is obtained, £ = 0.24 ±0.25. It may be argued that these surveys tend to 
select low-luminosity objects and the comparison with the SGP data could not 
be entirely significant, but an analysis on restricted absolute magnitude slices of 
the SGP sample shows no correlation of the clustering with the QSO absolute 
luminosity. Besides, the data of the Palomar Transit Grism Survey (Kundic 
1997, Stephens et al. 1997) allow measuring the amplitude of the TPCF at 
redshifts higher than 2.7 and the result, r Q = (18 ± 8)/i _1 Mpc, suggests that 
the trend of increasing clustering persists. 

If we describe the evolving correlation function in a standard way: £(r, z) = 
(r/r r 7 (l + z)-^ +t \ where e is an arbitrary (and not very physical) fitting 
parameter, we obtain e = —2.5 ± 1.0, which is inconsistent with the e ~ 0.8 
observed for faint galaxies at lower redshifts (Le Fevre et al. 1996, Carlberg et 
al. 1997, Villumsen et al. 1997). Great care should be exercised however when 
carrying out this comparison. Are the faint lower-redshift galaxies representative 
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of the same population of galaxies for which recent observations by Steidel et al 
(1998) show substantial clustering at z ~ 3.1 ? Are the Lyman-break galaxies 
progenitors of massive galaxies at the present epoch or precursors of present 
day cluster galaxies (Governato et al. 1998)? Can a mass scale of the order 
10 12 Mq be made consistent with the observations of field galaxy clustering at all 
redshifts (Moscardini et al. 1998) via an interplay between the clustering of mass 
(decreasing with increasing redshift) and the bias (increasing with redshift)? 

If we treat clustering according to the models in which quasars at z < 3 are 
associated with interactions, then we may relate their clustering properties - fol- 
lowing Matarrese et al. (1997) - to those of transient objects, which is definitely 
different from the case of galaxies which, depending on the physical scenario, 
can be associated with a merging or object-conserving paradigm of long-lived 
objects. In this way the observed clustering is the result of the convolution of 
the true clustering of the mass with the bias and redshift distribution of the 
objects and differences in each of these factors may lead to different results. If 
we think of QSOs as objects sparsely sampling halos with M > M m i n an increase 
of their clustering is expected because they are sampling rarer and rarer over- 
densities with increasing redshift. We may also ask what are the typical masses 
which allow reproducing the observed clustering. As we can see from Fig. || an 
-^min = 10 12 — 10 13 M & would provide the desired amount of clustering and 
evolution, and would be also in pleasing agreement with the predictions derived 
in Sect. 3 from the analysis of the QSO luminosity function. 
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